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ABSTRACT 


After  reviewing  some  of  the  previous  work  concerning 
kinetic  equations  for  desorption  a  comparison  is  made  bet¬ 
ween  the  Kramers-Fokker-Planck  equation,  derived  from  clas¬ 
sical  mechanics,  and  a  Fokker-Planck  type  equation  (GKTT 
equation)  recently  developed  by  Gortel,  Kreuzer  et  al,  from 
a  quantum  statistical  theory  of  phonon-mediated  desorption. 
To  facilitate  the  comparison  both  equations  are  expressed 
in  terms  of  scaled  action-angle  variables.  It  is  found  that 
by  ignoring  terms  of  the  order  of  the  thermal  fluctuations 
in  the  GKTT  equation,  it  reduces  to  the  low  friction  case 
of  Kramers'  equation.  The  friction  coefficient  is  identi¬ 
fied  and  when  expressed  in  terms  of  the  microscopic  para¬ 
meters  is  identical,  up  to  a  multiplicative  constant,  to 
the  expression  obtained  classically  by  Caroli  et  al. 
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I. 


INTRODUCTION 


A  gas  particle  is  said  to  physisorb  onto  the  surface 

of  a  solid  if  the  net  interaction  between  the  gas  particle 

and  the  solid,  accounted  for  by  an  effective  wall  or 

surface  potential  V  (r),  localizes  the  gas  particle  near 

s 

the  solid  surface.  The  surface  potential  typically 
consists  of  a  long-range  attractive  tail,  which  is  essen¬ 
tially  the  interaction  energy  between  the  mutually  induced 
fluctuating  dipole  moments  on  the  adsorbing  gas  particle 
as  well  as  in  the  solid,  and,  a  short-range  repulsive  core 
which  is  largely  due  to  increasing  charge  density  fluctua- 
tions  as  the  adsorbing  particle  becomes  confined  close  to 
the  surface.  Those  gas  particles  trapped  into  the  bound 
states  of  the  surface  potential  constitute  the  adsorbate. 
Typical  binding  energies  for  physisorption ,  as  measured  by 
the  heat  of  adsorption  Q,  are  less  than  or  even  much  less 
than  1  eV;  e.g.  for  Xe-W,  Q  =0.402  eV  and  for  He-C,  Q  = 
0.012  eV. 

Physisorption  is  one  end  of  a  range  of  adsorption 
phenomena  of  which  the  other  extreme  is  chemisorption. 

For  the  latter  the  adsorbed  molecule  undergoes  a  chemical 
change  by,  for  instance,  forming  bonding  orbitals  with  the 
solid.  For  chemisorption,  Q  is  typically  a  few  electron- 
volts.  In  this  thesis  we  shall  deal  only  with  physi¬ 
sorption  . 
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Theories  of  physisorption  fall  into  two  categories; 
equilibrium  and  nonequilibrium.  Whereas  equilibrium 
theories  attempt  to  calculate  quantities  such  as  the 
adsorption  isotherm,  specific  heats,  isosteric  heat  of 
adsorption  etc.,  nonequilibrium  theories  address  themselves 
to  the  kinetics  of  the  adsorption  and  desorption  processes. 
The  first  such  kinetic  theory  was  developed  at  a  phenomen¬ 
ological  level  by  Langmuir  in  1918'*'.  The  first  quantum 

mechanical  theory  was  put  forward  by  Lennard- Jones , 
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Devonshire  and  Strachan  around  1935  .  They  argued  that 

the  desorbing  gas  particle  acquires  the  energy  necessary 

to  leave  the  potential  trap  through  its  coupling  to  "the 

phonon  bath  of  the  solid.  This  point  of  view  has  been 

successfully  exploited  in  the  recent  quantum  mechanical 

3  4 

theory  of  Gortel,  Kreuzer  et  al  '  ,  where  second  order 

perturbation  theory  was  used  to  calculate  the  phonon- 
mediated  desorption  rates.  The  potentials  considered  in 
the  above  theory  developed,  at  most,  only  a  few  bound 
s  ta  te  s  . 

Potentials  with  a  large  number  of  states  and  class¬ 
ical  theories  have  also  been  put  forward.  For  gas-solid 
systems  at  low  coverage  one  argues  that  a  stochastic 
description  based  on  a  master  equation  is  sufficient.  Such 
master  equations  for  desorption  kinetics  have  been  cons¬ 
tructed  for  classical  models  by,  e.g.  Pagni"*,  Pagni  and 
6  7 

Keck  ,  Muller  and  Brenig  ,  etc.  A  theory  based  on  quantum 
statistical  mechanics  has  recently  been  developed  by 
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Kreuzer  et  al  '  '  .  Once  confronted  with  a  master  equa¬ 

tion  one  typically  works  at  deriving  a  simpler  partial 
differential  equation  of  the  Fokker-Planck  type.  This  has 
indeed  been  done  by  the  latter  authors,  thus  making  a 
connection  to  the  pioneering  work  done  by  Kramers  in  1940^. 

Although  the  papers  by  Kreuzer  et  al  derive  a  Fokker- 
Planck  type  equation,  which  they  use  to  study  desorption 
kinetics,  an  explicit  comparison  with  Kramers' work  was  not 
given.  This  is  the  task  of  this  thesis.  In  chapter  II  we 
will  briefly  review  the  Kramers  theory  and  later  develop¬ 
ments,  and  give  an  outline  of  the  quantum  statistical  theory 
of  Kreuzer  et  al .  To  facilitate  an  explicit  comparison  of 
these  two  theories,  we  will,  in  section  III. 1-3,  cast  both 
in  terms  of  action-angle  variables.  In  this  form  the 
quantum  statistical  theory  becomes  somewhat  more  trans¬ 
parent,  allowing  a  justification  of  some  of  the  earlier 
approximations,  and  indeed  to  go  beyond  them.  To  our  great 
surprise  we  found  that,  in  the  low  friction  limit,  the 
microscopic  expression  for  the  friction  coefficient  in  the 
Kramers-Fokker-Planck  equation,  as  derived  from  the  quantum 
statistical  theory,  is  identical,  up  to  a  multiplicative 

constant,  with  that  obtained  from  a  classical  theory  by 

12 

Caroli  et  al  .  This  is  explained  by  the  fact  that  low 
friction  results  from  weak  coupling  in  the  gas-solid  system 
and  hence  many  one-phonon  processes  must  occur  before  a  par¬ 
ticle  can  desorb,  a  situation  that  can  also  be  described  by 
classical  mechanics. 


II.  REVIEW  OF  PREVIOUS  WORK  ON  KINETIC  EQUATIONS  FOR 


DESORPTION 

II . 1  Classical  Theories 

In  1940  Kramers  developed  a  theory  of  a  particle 
trapped  in  a  potential  well  out  of  which  it  can  escape 
through  the  action  of  its  Brownian  motion^.  Intended  as 
a  microscopic  theory  of  chemical  reactions  it  has  found 
application  in  the  problem  of  gas  particles  desorbing  from 
a  solid  surface.  Kramers'  theory,  which  we  want  to  briefly 
outline  in  this  section,  requires  the  construction  of  a 
Fokker-Planck  type  equation  for  the  particle  distribution 
function  in  phase  space. 
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We  start  with  the  Kramers-Langevm  equation  , 


dp 

dt 


F(r) 


up  +  3  (t) 


(2.1) 


where  p  = mr  is  the  momentum  of  the  particle  of  mass  m 
trapped  in  the  potential,  r  its  coordinate,  F(r)  =- 


3r 

is  the  force  exerted  on  the  particle  by  the  potential,  and 
the  friction  coefficient  ri  is  related  to  the  random  force 
3  ( t)  by. 


n 


with 


1 

rnkBT 


<3(0)3  ( t)>  dt 

■ 

0 


<  3(0)  3  (t)  >  =  C6  (t) 


(2.2) 


(2.3) 


being  the  time  correlation  function  for  the  stochastic 
variable  3(t) .  From  this  stochastic  equation  an  equivalent 
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w 


Fokker-Planck  equation  may  be  derived  for  the  probability 
density  n(r,p,t).  If  the  particle  undergoes  a  random  walk 
in  phase  space,  then 


5 


n (r,p, t)  = 


g(r,p, t;r' ,p' ,t* )n(r' ,p' ,t' )  dr'  dp'  ,  (2.4) 


where  the  propagator  g ( r , p , t;  r '  ,  p '  ,  t 1 )  satisfies  the 
Smoluchowski -Chapman -Kolmogorov  equation 


g  ( r ,  p ,  t ;  r*  ,p'  ,t')  = 


g  (r  ,p,  t;r",p",t")  g(r",p",  t";r'  ,p'  ,t') 


dr"  dp"  . 


(2.5) 


To  derive  a  second  order  partial  differential  equa- 
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tion  we  follow  Kac  and  Logan  and  define  jump  moments  for 
small  time  increments  At  =  t-t"  : 


j  (r-r"  )  g  (r,p,  t;r",p",  t~At)drdp=  <(r-r")  >^_ 


-  Ar(r",p")  At,  (2.6a) 


(p-p")  g(r,p,  t;r",p",t-At)  drdp=  <(p-p")  > 


At 


A  (r",p")  At,  (2.6b) 


* 

(r-r")  (p-p")g(r,p,t;r",p"/t-At)drdp=<(r-rJ')(p-p")>At~B^(r"  ,p")  At  , 

(2.6c) 

(p-p")2g(r,p,t;r",p",t-At)drdp=<  (p-p")  2>  ~B  (r",p")At  ,  (2.6d) 

J  Xr  ir 

J  (r-r")2g(r,p,t;r",p",t-At)  drdp=<(r-r")  2>  ~  B^r  ( r  "  ,  p"  )  A  t  .  (2.6e) 

Now,  multiply  (2.5)  by  a  well-behaved  test  function,  cj)(r,p), 
and  integrate  over  variables  r  and  p  to  obtain. 
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g(r,p,  t;  r '  ,p'  ,  t)  (J)  (r,p)  dr  dp  = 


g(r,p,  t;  r"  ,p"  ,  t") 


g(r"  ,p"  ,  t";  r*  ,p'  ,  t* )  <f>  (r, p)  dr"  dp"  dr  dp 


(2.7) 


Expanding  <J)(r,p)  in  a  Taylor  series  about  (r",p") 
*tr,p)  =  <t>  (r"  ,p")  +  (r-r")  f|  +  <p-p")  || 


1  .  „,2  324> 

r=rU  2 !  1  3p3p 

p=p" 


r=r 

P=P' 


^  2  , 

+  (r-r") (p-p")  9 


3r  3p 


r=r 

p=p' 


1  /  ii  \ 2  31 

2  I  (p  P  }  3p3p 


r=r 

p=p" 


+  •  •  •  / 


(2.8) 


and  substituting  into  the  right  hand  side  of  (2.7)  gives. 


g(r,p,  t;  r*  ,p'  ,t)  (|)(r,p)  dr  dp 


g (r",p",t";  r'  ,p'  ,t')  <J)(r",p") dr"^" 


g  ( r " , p" ,  t" ;  r ' , p  Vt ' )  fer  [Ar  ( r ",p" )  A t  +  0  ( ( A t) 2)  ] 

2 

+  ^rr  [Ap(r",p'')At+0((At)2)]  +  -gp;f^r  [B^  (r”,p")  *t+ 0(  (At)  2)  ] 
2  2 

+  i-^4-  m  (r",p")At+0(( [B  (r\p")At+0((At)2)] 


2  3r»2  rr 


2  3p"2  PP 


+  0(  (At)  )  \  dr"  dp" 


(2.9) 


We  integrate  by  parts,  rearrange  indices,  and  divide 


by  At  to  obtain 


7 


_1_ 

At 


<Mr,p)  [g(r,p, t+At;r* ,p' , t' )  -g(r,p,t;r' ,p' ,t') ]dr  dp 


<f>  (r,p) 


9r  r 


[A  (r,p)  g(r  ,p,  t;r' ,p' ,t')  3  - 


9  A 

fAp(r'p)  g(r'P/t;r'/pVt')]  +  3^-  [B  (r,p)g(r,p,t;r*  ,p'  ,t')] 

92  1  a 

+ - 2  tBrT-(r'P)<?(r'P/t'r,'P,'t'  )  j  (r,p)g(r,p,t;r',p,,t’)] 

23r  J--L  +  dp^  pp 


+  O(At)  >■  dr  dp 


(2.10) 


Multiplying  both  sides  by  nCr'^p'^t'),  integrating  over  r* 
arid  p',  and  finally  letting  At 0  results  in 


*(r,p)  [^£LtL]dr  dp  = 


(f)  (r,p)  [Ar  (r,p)  n(r,p,  t)  ]  - 


[A  (r,p,t)n(r,p,t)  ]  +  j  [...]|drdp 


(2.11) 


and  so,  since  the  test  function  is  arbitrary 
9n  _ 

9r 

,2 


2  2 

[An]  -  7^-  [A  n]  +  i I -X?  [B  n]  +^~  [B  n] 
9t  9r  r  9p  p  2[g  2  rr  9r9P  rP 


+ 


9P 


o-  [B  n] 

2  pp 


(2.12) 


The  jump  moments  are  calculated  from  the  equations 
of  motion.  Integrating  p =  mr  gives 


t+A  t 


t+A  t 


m(r-r")  =m 


r  dt'= 


t 


p(t')  dt\ 


(2.13) 


and  so 


t+A  t 


<r-r">  =  — 
m 


<  P  (  t  '  )  >dt  '  S! 


£  At  , 
m 


(2.14) 
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which  gives 

A  ( r  "  ,  p  " )  ~  H 

r  m 

To  obtain  <(r-r")2>  we  square  (2.13) 


t+At  t+At 

p ( s) p ( s ' )  ds  ds  '  > 

=  -^2  <  (pAt)  (pAt)  > 
m 


m 


=  0  (  (At)  2 )  . 


< (r-r")  >  =■ 


< 


m 


(2.15) 


B 

rr 


(r" 


P")  ~  0 


(2.16) 


Similarly  we  may  integrate  (2.6a)  and  calculate  the  jump 
moments  involving  p,  although  it  is  necessary  to  use  (2.2) 


t+At 

and  the  fact  that  3 ( t)  is  random,  i.e.,  J 

t 

Substituting  these  moments  into  (2.12) 
the  Kramers  equation. 


9n  _  8 V(r)  8n  p  5n 


™  tr  +  h  [pn  +inkDT  r 

o t  3 io  3p  m  3 r  3p  B  3p 


&(t' )  dt'  =0. 
finally  gives 


(2.17) 


where  n  En(r,p,t)  . 

We  next  want  to  use  the  Kramers  equation  to  study 
physisorption  kinetics.  Because  a  gas  particle  physisorbed 
onto  a  solid  surface  is  rather  weakly  coupled  to  the  thermal 
excitations  of  the  solid,  it  will  perform  many  oscillations 
in  the  surface  potential  well  before  it  acquires  enough 
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energy  to  desorb. 


cient  r)  in  (2.17) 
that 


r\  <<  U) 


V 

u 


This  implies 
is  small,  i.e 


/ 


where 


0)  =  Y 

o 


that  the  friction  coeffi- 

20 

according  to  Suhl  et  al 

(2.18) 


(2.19) 


is  the  angular  frequency  of  oscillation  at  the  bottom  of  the 
well  of  depth  Uq  and  of  range  y  \  In  the  low  friction  limit  the 
particle  orbits  in  phase  space,  for  many  successive  oscill¬ 
ations,  are  within  the  thermal  fluctuations  in  r  and  p  and 
can  therefore  be  assumed  closed  for  all  practical  purposes. 
This  allows  us  to  introduce,  in  a  meaningful  way,  abtion 
and  angle  variables: 


p  dr 


(2.20) 


w  =  vt 


(2.21) 


Mathematically  speaking,  Kramers'  equation  cannot  be  written 
in  terms  of  the  action  and  angle  without  specifying  the 
potential.  However,  Kramers  has  argued  that  by  averaging 
(2.17)  over  an  orbit  in  phase  space  it  may  formally  be 
written  as 

ire  =  n  al7  [J"+  “V  v  Iy1  '  (2-22) 


where 

n  ( J,  t) 


o  n(J,w,t)  dw  . 


(2.23) 


Let  us  emphasize  that  (2.23)  is  only  a  formal  relation. 
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To  study  (2.22)  further,  in  the  context  of  physisorption 
kinetics,  we  need  the  relation  between  the  action  variable 
J  and  the  frequency  v.  It  thus  seems  worthwhile  to  connect 
the  action  variable  J  to  the  energy  E  of  the  particle 
trapped  in  the  surface  potential.  This  can  only  be  done 

after  a  potential  V(r)  has  been  specified. 

15  n 

Bak  and  Andersen  choose  a  potential  V(r)  =c|r| 

and  show  that  E  aJv.  Eliminating  the  action  in  favour  of 

the  energy  in  (2.22)  and  choosing  boundary  conditions 

appropriate  for  desorption,  as  well  as  an  ansatz  n(E,t)  = 

n(E)  exp  (-At) ,  they  have  shown  that  n(E)  can  be  expressed 

in  terms  of  confluent  hypergeometric  functions  and  that  A  ^ 

can  be  identified  as  the  desorption  time.  We  will  present 

such  an  analysis  in  chapter  III  for  the  more  realistic 

Morse  potential, 

-2y(r-r  )  -y (r-r  ) 

V  ( r)  =  U  (e  °  -  2  e  °),  (2.24) 

o 


after  we  have  derived  (2.22)  from  the  quantum  statistical 
theory . 

12 

Recently,  Caroli  et  al  have  derived  the  Kramers- 
Langevin  equation  (2.1)  starting  from  a  microscopic  model 
based  on  classical  mechanics.  This  theory  yields,  for  a 
Morse  potential,  a  microscopic  expression  for  the  friction 
coefficient;  namely 


4  a) 


m 


D  M 


u  2 
(— ) 


n 


r 


/ 


(2.25) 


2 

where  u^  =  U^/ftco^  ,  r  =  2moo^/hy  ,  go  is  the  Debye  frequency 

of  the  solid,  and  M  is  the  mass  of  a  solid  atom. 

s 


II •  2  Quantum  Statistical  Theory 

We  will  now  outline  the  quantum  statistical  theory 
of  physisorption  as  developed  by  Kreuzer  et  al^'^'^. 
According  to  quantum  mechanics  the  surface  potential  will 
develop  a  number  of  bound  states  at  energies  which  may 
be  occupied  by  an  adsorbed  gas  particle.  The  probability 
o£  the  gas  particle  being  in  energy  state  is  given  by 
the  occupation  function  n^.  It  is  assumed  that  the  phonon- 
mediated  transitions  between  different  bound  states  as  well 
as  to  and  from  the  gas  particle  continuum  may  be  regarded 
as  a  Markovian  process  and  hence  the  occupation  functions 
n^  are,  at  low  coverage,  subject  to  a  set  of  rate  equations 

dn. 

1 

dt~ 


l  R .  .  n  .  -  l 

k  ±3  3  jk 


R.  .  n. 
Ji  i 


(2.26) 


where  the  indices  i  and  j  run  over  all  bound  states  and 

continuum  states  accessible  to  the  gas  particle.  The 

transition  probabilities  R^  have  been  calculated  by  Gortel 

8 

Kreuzer  and  Teshima  for  a  one  dimensional  Morse  potential 
(2.24)  using  second  order  perturbation  theory. 

For  gas-solid  systems  such  as  Xe-W  there  are  over  a 
hundred  bound  states.  It  therefore  seems  reasonable  to 
replace  the  sums  over  i  and  j  by  integrals  over  the  scaled 


energy 


12 


£  =  E/iiu)  .  (2.27) 

Letting  the  occupation  function  go  over  into  a  function  of 
energy , 


ni(t)  -*n  (e±/ 1)  ,  (2.28) 

9 

we  may  now  write  (2.26)  as  a  continuous  master  equation  , 
3n(£,t) 


3 1 


de*  p(e')W(E,£,)n(e',t)  - 


. 

-u 


-u 


de'  p  (  e  ' )  W  (  e  '  ,  e)  n  ( £  ,  t)  , 


(2.29) 


where 


P(e)  =  — —  .  (2.30) 

/yir 

V 

For  the  bound  state-bound  state  transitions,  i.e.  for 

9 

-u  <  e , e 1  <  0,  Kreuzer  et  al  give  the  integration  kernel 
as 


W  (  e  ,  e  1 )  =  3tt/t  m  ^/ee  '  (  e-e  1  )  ^  [exp  [6  (  e-e  1  )  J  ] 


-1 


X 


0 (1-E+E ' ) 0 ( £“£ ’ ) 


T(/ru  - /-r.£  +  ~-)r(/ru  + /-rs  +  _•) 

O  O  2  O  2 

_ 


T(/ru  - /-r£  1  +T)r(/ru  + /-r£  +  ~-) 
o  2  o  2 


+  0  ( 1“ £ ' +E)  0  (  £ ' — £ ) 


r  ( /ru  -  /-r£  '  +7r)r(/ru  + /-r£  '  +  tt) 

O  2  O  2 


T(/ru  - /-r£  +  ^-)T(/ru  + /-r£  +  ^-) 

O  2  O  2 


(2.31) 


where  0  is  the  step  function  and  6  =  'fta^/k^T. 

They  also  give  the  integration  kernel  for  bound  state- 

continuum  transitions,  i.e.  for  -u  <  e'  <  0  and  0  <  e  <  °° 

o 
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3tt 


W(e,e')  =  /r  oo 

2  DM 


m 


/-££ '  0 (1+e') 0 (l+e'-e) 


X  (exp[S  <£-£')]  -l)_1  (£-£')  - sinh  (2tt/tF) 


sinh  [ it /re ]  +  cos  [u/ru  ] 

o 


x 


r  (./ru  +i/rF  +  =-) 
O  2 


r  ( /HT  +  /-re  '  +  j)  r  ( /ruQ  -  /-re  1  +  j) 


(2.32) 


and  for  continuum-bound  state  transitions 


W(e',e)  =  T  ^  WD  if  >/zee~r  0(l+e')  0(l+e'-e) 


x| (exp [6 (e-e ' ) ] -l)  1+l|(e-e') 


sinh  (27r/re) 

2  ,  _  . — ,  2 


T  (  /ru  +  i  /re  + 
o  2 


T  (/ru  +  /-re  '  +  4)  T  ( /ru  -/-re1  +  4) 
o  Z  o  2 


sinh  (tt/tFj+cos  (tt^tu^) 


(2.33) 


The  integration  kernel  W(e,e')  is  plotted  in  figures 
1  and  2.  Examining  these  figures  we  see  that  the  bound 
state-bound  state  transition  rates  are  much  larger  than  the 
bound  state-continuum  transitions  for  most  energy  values. 
Physically  this  means  that  an  adparticle  will  undergo  many 
bound  state-bound  state  transitions  before  desorbing.  This 
large  separation  in  time  scales  leads  us  to  believe  that  the 
adparticle  is  executing  a  random  walk  amongst  the  bound  state 
energy  levels.  We  suspect  that  a  Fokker-Planck  equation  will 
describe  the  adparticle  kinetics.  We  also  note  that  near  the 
top  of  the  potential  (e=0)  the  bound  state-bound  state  and 
bound  state  continuum  transition  rates  are  of  the  same  order 
of  magnitude.  We  do  not  expect  a  Fokker-Planck  type  equation 
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Fig.  la  and  Fig.  lb  . 


Hie  integration  kernel 
W  ( e , £ ' )  for  the  full 
range  of  £  and  e1.  Transitions  Into 
Courtesy  of  Prof.  Higher 

Kreuzer. 


Bound  States 


-1.0 
-23.0 


Transitions 

Into 

Continuum 


Transitions  Into 
Lower 

Bound  States 


-1.0 


-23.0 


1.0 


Transitions 

From 

Continuum 


Ficr.  lb 
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to  be  valid  in  this  region. 


Converting  the  integral  operators  in  (2.28)  into 
differential  operators  will  indeed  yield  a  Fokker-Planck 
equation.  Toward  this  end  we  define  a  propagator  g(e,t;e't' ) 

m 

by 


n(e,t)  = 


de'  g(e,  t;  e  '  ,  t' )  n  (e  '  ,  t' )  , 


(2.34) 


-u 


which,  for  small  time  increments  At  =  t-t' ,  can  be  expanded 


as 


g(e,  t+At;  e  '  ,  t)  =  6(e-e')  +  At-f  p  (e)  W(e  ,  e  ' )  -S(e-e') 


de" 


-u 


o 


p  (  £  "  )  W  (  £  "  ,  £  1  ) 


(2.35) 


It  also  satisfies  the  Smoluchowski-Chapman-Kolmogorov 
equation 


g(e,  t;  e  1  ,t' ) 


de"  g  ( e  ,  t;  e"  ,t")  g(  e" ,  t"  ;  e'^t')  .  (2.36) 


Following  the  procedure  used  to  develop  the  Kramers  equa¬ 
tion,  we  multiply  both  sides  of  (2.36)  by  a  well-behaved 
test  function  t(e)  and  integrate  with  respect  to  e ,  to 
obtain , 


t  (e)  p  (e)  g(e,t;e,,t')  de  = 


-u 


CO 

r 


-U 


p  (e)  t  (e)  g  (e  ,  t;e",t")  g(e"  ,  t"  ;e',t') 
,  -de  de"  (2.37) 


The  test  function  may  be  written  as 


T(e)  =  l  ■—  (e-£")n  x(n)  (£")  , 


n=0 


n 


(2.38) 


which,  upon  substitution  into  the  right  hand  side  of  (2.37) 
gives , 
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-u 


x(e)  p  ( e)  g  (  e  ,  t;e'  ,t' )  de  = 

-u 


+  At  [  ± 

n=l 


n! 


x(e)  p(E)g(E,t";E'  ,t' )  de  + 


(n) 


an  (e  ")  g(e  ,  t;e',t' )  t  (  e")  de"  , 


(2.39) 


where,  for  small  time  intervals  At  =  t-t", 


“n(e">  =/t 


n 


dep(e)(e-e")  g  ( e ,  t+At;  e  ",  t) 


p(£,,)  j  P  (£)  (e~e")n  w(e,e")de=p  (en)an(e")  .  (2.40) 


o 

Rearranging  terms  and  performing  some  necessary  partial 

» 

integration  in  (2.39)  leaves  us  with. 


t  (e)  p  (e) 


g(e,t;e',tl)  -  g  ( e  ,  t-At;  e  1  ,  t 1  ) 


At 


de 


•  r  vAV 

T(£)p(£)|  l 

^n=l 


°°  n  n 

(-1)  8 


I  _  (a  (e) g(e, t?e ' ,t' ) )  >  de  . 
n!  9£n  n  / 


(2.41) 


Taking  the  limit  At->0  and  multiplying  (2.41)  by  n(e',t') 
and  integrating  over  e'  gives,  using  (2.34), 

T^e)|Ft  "  l  - -  (an(e)n(e,t)  ) 

Since  (2.42)  holds  for  any  test  function  x(e),  the  expres¬ 
sion  in  brackets  is  identically  zero.  Hence,  (2.28)  may  be 
rewritten  as  a  Kramers-Moyal  expansion, 

°o  n  n 

[  p(  e)  n  (  e,  t)  ]  =  l  - [a  (e)  p  (e)  n  (e  ,  t)  ]  .  (2.43) 

9t  n=1  n  I  9en  n 


\de  =  0  . 

J  (2.42) 
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Truncating  the  series  at  n=2  gives  a  Fokker-Planck  type 
equation , 


g 

[p  (e)  n  (  e  ,  t)  ] 


_9_ 

9  e 


( ot1  (e)  p (e) n (e,  t)  ) 


+ 


+  \  ~^~2  (<*2  (e)  p(e)n(e,t)  )  *  (2.44) 

9  e 

We  shall  refer  to  this  equation,  with  or  without  approxima¬ 
tions  to  the  moments  ot^(e)  and  a^{e) ,  as  the  GKTT  equation. 

Kreuzer  et  al  found  approximate,  analytic  expres¬ 
sions  for  the  moments  and  in  the  temperature  regime. 


5  _1_ 
3  u 

o 


<  < 


6 


<  < 


(2.45) 


namely, 

a1(e)  *  +72ti  a  /+re  a'5  -  1]  ,  (2.46) 

s 


a2U)  s  6,1  “d  m 


m  , -  “5  r24  120  2520B. 

/-re  A  [~r  + — o - - 5 — J  / 


5  '  A2 


A' 


(2.47) 


where 


u  /  -e 
o 


0  tanh  ^  G 


u 


(2.48a) 


B  =  J_  /Z  (  UQ)j  1  - 

4uo/uo  e  h+e 


rU 


0  tanh  /— 


-e 


u. 


u 


(2.48b) 


Using  these  definitions  and  (2.44),  Kreuzer  et  al  were  able 
to  derive  an  approximate  expression  for  the  desorption 


..  10 

time  , 


,  M  2  _  6u 

1  -1  s  r  1  o 

t ,  ~  -  a)  —  —  -jz -  e  , 

d  72tt  D  m  „2  6u 

uo  0 


(2.49) 
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and  a  macroscopic  law 


9 


at 


e 


ot  (e) 


(2.50) 


where 


0 


£  = 


e  n  ( e  ,  t)  p  ( £)  d£ 


-u 


is  the  average  energy  of  the  adparticle 


(2.51) 


III.  the  gktt  equation  in  action-angle  variables 


III.l  Action-Angle  Variables  for  a  Morse  Potential 

The  major  task  of  this  thesis  is  to  relate  the 
Kramers  equation  (2.22)  used  in  the  classical  theories  of 
desorption,  to  the  GKTT  equation  (2.44),  which  was  derived 
from  a  quantum  mechanical  theory  of  physisorption .  To 
facilitate  such  a  comparison  we  first  introduce  into  (2.44) 
action  and  angle  variables.  We  recall  that  this  equation 
was  derived  for  gas  particles  trapped  in  a  one  dimensional 
Morse  potential  (2.24),  which  we  prefer  now  to  write  as 


-y (r-r  ) 


V  (r)  =  U  (e 
o 


by  adding  a  constant  UQ  to 


Hamilton- Jacobi  equation 
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in  the  potential  (3.1) 

(1^)  +  V (r)  =  E 
2m  dr 


2 

-  1)  (3.1) 

(2.24).  Starting  from  the 
for  a  gas  particle  of  mass  m 


(3.2) 


where  W=W(r,p)  is  Hamilton's  characteristic  function,  we 
introduce  in  the  standard  way,  the  canonical  momentum. 


P 


3W 

3r 


(3.3) 


This  may  be  inserted  into  the  definition  for  the  action 
variable  (2.20)  to  obtain 

3  W 


j  =  o  p  dr  =  o 


3r 


dr 


, _  /  UQ  -Y(r-r0)  2 

=  /2mE  4>  J 1  ~~  (e  -  1)  dr 

hj 


(3.4) 
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I  1 1 ) 
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-y(r-r  ) 

Setting  y  =e  we  get. 


J  =  - 


/2mE 


Y 


U 


U. 


°  1  /Uo  -|  \  1 

r+2  E“  y  -  (X_1)  ^2  dy 


C)  - 


2  7T 


„  /2m  [/u~  -  /u~  -E  J 
Y  o  o 


The  action  ranges  in  value  from  J  =  0  at  the  bottom  of 
the  well  (E  =0)  to 

2tt 


(3.5) 


J  =  J  =  —  /2mU 
o  y  o 


(3.6) 


at  the  top,  where  E  =U  .  In  dimensionless  units  we  have 

o 


J  =  h 
o 


/2mU 
_ o_ 

*hY 


=  ha. 


(3.7) 


We  may  invert  equation  (3.5)  to  express  E  as  a  function 
of  J 


E  _  Y  j  _  (_Y_)2 

tt  J  2m  J  2tt  2m 


(3.8) 


A  word  of  warning  is  in  order  here.  The  integral  in  (3.4) 
is  not  defined  when  E  =  U  ,  i.e.  when  the  particle  is 
no  longer  confined  to  the  surface  region.  It  will  be  seen 
later  that  this  introduces  singularities  into  the  Fokker- 
Planck  (GKTT)  equation.  However,  the  singularity  structure 
of  the  GKTT  equation  with  suitably  chosen  expressions  for 
the  moments  doesn't  change. 

We  can  now  find  the  frequency  of  oscillation  in  the 
Morse  potential  by  differentiating  the  energy  with  respect 
to  the  action. 
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v  = 


dE 

dJ 


y/u~ 

o 


tt/2 


m 


-  (^-)  ^ 

v  2tt;  m 


(3.9) 


Using  (3.8)  the  action  J  may  be  eliminated  in  favour  of 
the  energy  E, 


v 


(3.10) 


Note  that  this  frequency  is  energy  dependent  as  it  must 
be  for  any  realistic  potential  that  develops  both  bound 
states  and  continuum  states.  At  the  bottom  of  the  well  the 
frequency  is 


(3.11) 


We  check  this  result  by  examining  the  small  oscillations 
of  the  particle  about  the  well  minimum.  Energy  conserva¬ 
tion  gives  us. 


E  = 

2m 


+  V(r)  =  E-  +  u  (e 
2m  o 


-y (r-rQ) 


-1) 


(3.12) 


For  small  displacements  (r-r  )  from  r^  we  may  write  the 
exponential  term  as 


-Y (r-r  ) 

e  z  1  -  y (r-rQ)  , 


so  that 


E 


2 

„  P_  ,  5L 
2m  2 


2 


(r-r 

o 


) 


(3.13) 


(3.14) 


This  is  the  energy  of  a  particle  executing  harmonic  motion 
with  frequency 


*o  TiiOVBi  nx  od  *fi*  t  noiJo*  (  .  »ni"  J 
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v  = 


0) 
2  7T 


_1_ 
2TT  yj 


2U  Y 
o 

m 


fu~ 

=  1  1-9 

it  \/  2m 


(3.11) 


matching  the  result  obtained  for  the  Morse  potential  at 
E  =0.  In  closing  this  section  let  us  note  that  the  angle 
variable  w,  conjugate  to  the  action  J,  is 


,  ^  1  /2y 

W  =  Vt  +  W  =  77—  / — - 

o  2i  J  m 


(U  -E)  +w  , 
o  o 


(3.15) 


where  w  is  a  constant, 
o 


III. 2  Transformation  of  the  GKTT  Equation  to  Action-Angle 
Variables  * 

Before  transforming  the  GKTT  equation  to  action- 
angle  variables  we  introduce  a  dimensionless  action 


x  = 


(3.16) 


where 


2tt 


J  =  —  /  2mU  =  ha 
o  y  o  c 


(3.7) 


is  the  action  variable  at  the  top  of  the  Morse  potential 
well.  The  frequency  then  is 


v  =  —  ^  ( 1 — x ) 

7i  r 


/  u 

D  /  o  x 

—  / —  ( 1-x)  , 
7T  v  r 


(3.17) 


and  the  energy  is  given  by 

2 


e  =  -  —  (1-x) ^  =  -  u  (1-x) 2 


(3.18) 
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For  future  reference  we  note  that  the  Kramers  equation 
(2.22)  in  terms  of  the  variable  x  reads 


3n 

3t 


n 


i  a  rx  an,  a  r  - , 
—  U — -  — J  +^—  [xn] 


26u  ax  1-x  3x  ax 
o 


(3.19) 


To  transform  the  GKTT  equation  (2.44)  we  first  differen¬ 
tiate  equation  (3.18)  with  respect  to  x  and  get 
3e 


15  =  2uo(1  -x)  ' 


(3.20) 


so  that  the  partial  derivatives  with  respect  to  s  transform 


as 


and 


3f  _  3x  3f_ 
3e  _  3e  8x 


af 


2u  (1-x)  ax 
o 


anf 


n 


a  e 


n 


[12  Aj  f 

Lae  axJ 


n 


=  [ 


2u  1-x  3x 
o 


]  f 


-L  )  (-1)  v  {-7\K  l  O 

I^(2u)nk=l  (k-UMn-k)!  ^ 


n  n 


( 2n-l-k)  !  a 


(2uo) 


(3.21a) 


[■ 


(1-x) 


2n-l-k 


] 


(3.21b) 


Using  (3.21a)  the  GKTT  equation  (2.44)  reads,  in  the 
variable  x, 

TE  [p  (x)n(x,t)  ]  =^2  l^x  9x  {l^x  8x  [a2  (x)  P  (x)  n  (x' ] 

o 


i  i  a 


2u  1-x  3x  1 
o 


[a.  (x)p(x)n(x,t)  ]  ,  (3.22a) 


or  since 
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P  (x)  = 


/u-  (1-x) 
o 


(3.23) 


9n (x, t) 
9 1 


8u 


2  9x  1-x  9x 


i  a  a  (x)n(x,t) 

3  t— -■< — ] 


(1-x) 


2u  9x 
o 


,  a,  (x)  n(x,t) 

3  I  —  Tt  TIT - ] 


(1-X) 


(3.22b) 


Alternatively,  we  may  write  (3.22)  as, 
9  f  1 


a  ot9(x)  -  ,  a  (x)f  a9  (x)  -  (1-x)  a9  (x) 

[— - -  Ml  ^ -  f] 


9t  o,2  9x  ,2  9x  9x  2u  (1-x) 

ouq  (1-x)  o 

an  (x)  9  (x) 

L2u  (1-x)  _  2  ,4J  r  ' 

o  8u  (1-x) 

o 

where 


2  3 

8u  (1-x) 
o 


(3.24) 


f  =  f(x,t)  =  p(x)n(x,t) 


(3.25) 


III. 3  Some  Comments  on  the  Kramers-Moyal  Expansion 

9 

The  GKTT  equation  was  derived  from  the  master 
equation  (2.28)  by  a  truncation  of  the  Kramers-Moyal 
expansion  (2.43)  at  n  =2.  We  want  to  examine  briefly  the 
question  whether  a  Kramers-Moyal  expansion  in  the  variable 
x  would  after  truncation  also  lead  to  equation  (3.22). 
Using  (3.21b)  we  get  from  (2.43), 


oo  n  r 

9  r/  \  /  j_\i  V  "I  1  J 

at  [p(x)n(x't)]  =  >  \ 

n=l 


(-1) 


n  n 


n 


l  (~2) 


(2uq)  k=l 


(2n-l-k)  I  9' 


an  (x)  p  (x)  n  (x,  t) 


(k-1)  !  (n-k)  I  8xk 


(1-x) 


2n-l-k 


(3.26) 
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Substituting  (3.23)  into  (3.26),  multiplying  by  (1-x)  and 
interchanging  the  n  and  k  summations  gives 


3n(x,t)  "  (-l)k  3k  ,  ,  ,  .  ,  , 

- 3t  ,  TT  TT  tek(x)n(x,t)] 

k — 1  o  x 


(3.27) 


where 


6k(x)  =k  2 


( 2n-l-k)  ! 


a  (x) 
n 


V  _ 

L ,  lo2n  ,  ...  n  n  .  2n-k 

n=k  n !  2  (n-k)  I  uo  (1-x) 


(3.28) 


This  is  a  power  series  in  —  which  suggests  that  for 

uo 

UQ  >>  1  a  truncation  at  n  = 2  is  again  possible,  justifying 
equations  (3.22)  and  (3.24) .  The  condition  uq  >>1,  or, 

Uq  >>-fuj o  ,  implies  that  the  desorbing  gas  particles  must 

t 

undergo  a  great  number  of  individual  one-phonon  absorption 
processes  before  it  leaves  the  surface  potential  well. 

Note  also  that  the  thermal  fluctuations  in  the  x 
variable,  calculated  with  a  Maxwell-Boltzmann  distribution, 
are 

2  2  1 

<x  >  -  <x>  =  =  -g-  .  (3.29) 

ou  U 

o  o 

With  the  inherent  restriction  k  T  <<tiio  ,  needed  to  justify 

13  D  # 

the  one-phonon  calculation  of  the  transition  probabilities 
W(e,e')  (2.31),  we  see  that  the  condition  uq  >>  1  implies 

that  the  fluctuations  should  be  small.  This  is  important 
to  keep  in  mind  at  the  stage  where  we  will  try  to  appro¬ 
ximate  the  moments  a  (x) . 

n 

It  is  not  known  what  relationship  this  expansion 
has  to  the  method  of  expanding  in  a  parameter  ft  ,  where 
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is  large,  used  by  van  Kampen^^'^^  to  justify  the  trunca¬ 
tion  of  a  Kramers-Moyal  expansion. 


III. 4  Approximation  of  the  Moments  a  (x) 

- — - n- — — 

Approximate,  analytic  expressions  for  (x)  and 
a2 (x)  •  limit  of  large  r,  have  been  obtained  in 

ref.  [9]  and  appear,  respectively,  as  (2.46)  and  (2.47). 
In  the  x  variable  these  equations  read 


a1  (x) 


m  uo  , ,  >6  r,  ,-1,,  .  .  J  r10B 

+  72tt  w  - 2  d“x)  [tanh  (1-x)  ]  [-rr- 

s  r  0/1 


-5 


-1],  (3.30) 


and 


a2  (x) 


~  6TT03 


m  uo  ,,  ,6  r,  ,-1,.  ,  r5  r24  1206 


DM 


2  d-x)  [tanh  (1-x)  j  [—  + — 2 


252  OB 


]  , 


where 


A 


I  r  tanh  ^ ( 1-x) 
/u  o  (1-x) 


(3.31) 

(3.32) 


and 


1  fr  1  r  1  _  tanh  (1-x)  , 

4u  /u  1-x  ,  ,  ,2  (1-x) 

o J  o  l-(l-x) 


(3.33) 


To  approximate  a^(x)  and  (x)  ,  for  intermediate  values  of 
x,  by  yet  simpler  expressions  we  resort  to  examining  some 
numerical  examples.  From  figures  3  and  4  we  see  that 
a^(x)  and  (x)  are  nearly  quadratic  in  x  away  from  x=0 
and  x=l.  They  are  both  off  centre  with  a  degree  of 
asymmetry  depending  on  6;  a^(x)  becomes  more  asymmetric 
for  decreasing  6  while  the  opposite  is  true  for  a  (x) . 

Both  of  the  moments  go  to  zero  at  x  = 1  and. become  very 
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small  at  x  — 0.  The  moment  a^(x)  crosses  the  x  axis  for  a 
small  positive  value  of  x. 


and 


We  propose  the  following  approximation, 

a1(x)  ~  -72tt  u)d  - - j  (x-a) (1-x)  , 

s  r 


m  "1 

a2  (x)  ~  144tt  03d  - - 2  5"  (x+b)  (1-x)  . 

s  r 


(3.34) 


(3.35) 


Note  that  the  zeros  a  and  -b  are  of  the  order  of  the  ther- 

1 


mal  fluctuations. 


2a  (x) 


Sa2 (x) 


6u0 

x-a 

x+b 


Examining  figure  5  shows  that 

(3.36) 


is  qualitatively  the  correct  form  of  the  ratio  of  the 
moments.  The  correctness  of  the  factor  (1-x)  can  be  seen 
by  examining  figures  3  and  4,  where  we  plot 


a(x)  =  a(x)p(x)  =  - - -  .  (3.37) 

/u~  (1-x) 
o 

If  a^(x)  and  (x)  did  not  contain  the  factor  (1-x)  then 
the  a-^(x)  and  a2  (x)  would  not  remain  finite  at  x=l.  Figures 
3  and  4  also  show  that  a^(x)  and  $2(x)  are  roughly  linear, 
with  zeros  near  x=0. 


III. 5  The  New  GKTT  Equation 

Substituting  (3.34)  and  (3.35)  into  (3.22b)  gives 


9n  (x,  t)  _  J) 
9t 


1  x+b  3n(x,t) 


3x(26uo  1-x 


3x 


+  ( x  -  a  + 


2Su  (1-x) 
o 


-)  n  (x,  t) 


(3.38) 
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Fig.  5.  A  plot  of  the  ratio  2a1 (x) /6a2 (x) .  The  curve 


V  = 


x  -  .02 


x 


(a=J.023,  b  =  0)  matches  the  curve 


for  6=2,  to  within  a  few  percent.  It  is  not  shown 
since  it  lies  on  the  curve  for  6=2  most  of  the 
time 
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where 


(3.39) 


This  is  the  new  version  of  the  GKTT  equation.  Since 


a  and  b  are  of  the  order  of  the  thermal  fluctuations  in  x, 
which  are  small  by  virtue  of  (2.45) ,  we  may  make  further 
approximations,  namely 


x  +  b  ~  x 


and 


1 


x  -  a  + 


26u  (1-x)  ~ 

o 


where  it  is  understood  that  x  does  not  approach  the  values 
0  or  1.  These  approximations  yield. 


9n  (x,  t) 
8 1 


(3.40) 


This  is  identical  to  Kramers'  classical  equation  (3.19). 

By  making  appropriate  approximations  to  the  quantum 
statistical  Fokker-Planck  equation  (2.44)  [GKTT  equation], 
we  have  shown  that  it  is  equivalent  to  Kramers'  equation 
(3.19)  provided  that  thermal  fluctuations  are  ignored. 

This  constitutes  a  derivation  of  a  classical  equation  from 
a  more  basic  quantum  statistical  theory. 


Now  that  the  GKTT  equation  has  been  related  to 


the  Kramers'  equation  we  may  formally  identify  r\  as  a 
friction  coefficient.  Before  proceeding  we  wish  to  examine 


an  alternate  form  of  (3.40).  Starting  with 
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3f 

at 


1  a  a2(x)  at  a  a2(x)f  a2  (x)  -  (1-x)  a2 (x) 

0.2  ax  ,2  axJ  ax  [2u  (i-x)  +  0  2 ,3  ] 

°uo  ( 1— x)  o  8uq(1-x) 


+  [- 


(x) 


2u  (1-x) 
o 


,  a  (x) 

+  ~~2  - - 4  ]  f 

8u^  (1-xp 


(3.24) 


where  f  is  defined  by  (3.25),  we  substitute  in  equations 
(3.34)  and  (3.35)  to  obtain,  after  ignoring  the  thermal 
fluctuations , 


3_f  _  f  1  _3_ 

at  \26u  ax 

^  o 


r  x  a  f  a  r  4-1 

[Kf] 


(3.41) 


By  using  (3.7),  (3.11),  (3.13)  and  (3.14)  we  write 

equation  (3.41)  in  terms  of  the  action  J  and  frequency 
of  oscillation  v. 


a f  _  /  a  rj.  „  a f ^  a  j  vo  _ 

at  n(aj  [vkBT  ajJ  aj  (J  }  jq  v  f 


(3.42) 


It  is  readily  seen  that  (3.42)  contains  no  quantities 
specific  to  the  Morse  potential  and  so 'should  hold  for 
any  realistic  surface  potential  appropriate  for  physi- 
sorption . 
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HI-6  The  Macroscopic  Law  of  Friction  and  the  Friction 
Coefficient 

Kreuzer  et  al  have  obtained  an  approximate  macros- 
9 

copic  law  from  (2.44),  namely, 

(2.50) 

where  e  is  the  average  scaled  energy.  This  can  be  re¬ 
written,  using  (3.31),  (3.15),  and  (3.26)  as  follows: 


L.H.S  .  = 


—  =—!#=  2u„ (1-x)  , 


^  ^  9  x  ^  ^ 


at 


(3.43) 


where 


u  2 

R.H.S.-  a,  (x)  =  -  7  2tt  a)  (-?-)  u  (1-x) (x-a) 

1  D  M  r  o  ‘ 


=  -  2u  n (1-x)  (x-a) 
o 


(3.44) 


and  so  equating  both  sides  and  cancelling  common  factors 
gives  us, 

|^r  =  -  n  (x-a)  .  (3.45) 

This  is  the  macroscopic  law  of  friction  with  friction 
coefficient  g .  The  solution  to  equation  (3.42)  is 

x  =  a  +  (xQ-a)  e  ^  ,  (3.46) 

where  xq  is  the  initial  action  of  the  particle.  Physically 
this  means  that  if  a  gas  particle  is  put  into  a  bound  state 
x=xq,  it  will  lose  or  gain  energy,  exponentially,  until 
its  average  value  of  x  is  x=a,  where  a  is  of  the  order  of 
a  fluctuation  from  the  bottom  of  the  surface  potential  well. 
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This  is  the  same  interpretation  put  forward  in  ref.  [9] 
Note  that  the  microscopic  expression  (3.39)  for 
the  friction  coefficient. 


r]  =  36tt  a) 


_  u„ 
m  o 


D  M 


s  r 


=  9tt 


2  4 
U  y 
o 

M  3 
m  M  wn 
s  D 


(3.39) 


is,  up  to  a  constant  the  same  as  that  obtained  by  Caroli 

et  al  in  their  microscopic  model  based  on  classical 
12 

mechanics  .  This  agreement  of  quantum  and  classical 

theories  seems  to  stem  from  the  condition  u  >>  1  for  the 
/  ° 

validity  of  the  equations  (3.22)  and  (3.24)  which  implies 
that  a  large  number  of  one-phonon  transitions  are  necessary 
to  desorb  a  particle  from  the  bottom  of  the  surface  poten¬ 
tial  well. 

We  can  now  calculate  the  friction  coefficient  n  for 

a  particular  gas-solid  system,  which  we  choose  to  be  Xe-W, 

to  see  if  the  low  friction  assumption  is  indeed  satisfied. 

m 


2  13  -1 

With  the  parameters  ooD  =  5 . 3  x  10  s  ,  M 


•=  .714, 


s 

r  =  4969 ,  and  u  =11.56,  we  find  ri  =  2.3  x  lO'*'8  s  To  check 

o 

whether  this  value  is  appropriate  for  the  low  friction 
limit  we  rewrite  n  in  terms  of  the  oscillation  frequency 
at  the  bottom  of  potential,  i.e.  oo  q  ,  to  get 


0  m  f%3/2  ,  kBT, 

n  =  187r  (T)  6uo  (%  lf-) 


(3.47) 


hgT 

For  Xe-W  this  implies  n  s:  0.1  (w  -£—)  and  hence  n  << 


o  u. 


o  n 


which  is  indeed  Kramers'  condition  for  low  friction 
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HI*7  Desorption  Kinetics 

In  this  final  section  we  shall  study  desorption 
kinetics  based  on 


3n  _  _9_  f  1  x+b  9n 

3t  ’’  dx  I26u_  1-x  8x 


+  (x  -  a  + 


26u  (1-x) 
o 


-)  n 


(3.38) 


Before  examining  the  time  dependent  solutions  of  (3.38)  we 
first  find  its  equilibrium  solution.  Setting  =  0  and 
ignoring  terms  of  the  order  of  the  thermal  fluctuations 
reduces  (3.38)  to 


3  [:  1 


3x  2Suq  1-x  dx 


x  dn  ,  n 

t: —  +  xn  ]  ~  0 


(3.48) 


The  solution  to  (3.48)  is 

6u  (1-x) 2  6u  (1-x) 

n(x)  s  c1  e  +  c2  e 


,  -6u  (1-x) 

1-x  o  , 

-  e  dx  , 

x 


(3.49) 


where  the  constants  are  determined  by  the  boundary  condi¬ 
tions.  Setting  c^  = 0  gives 

6u  (1-x) 2  _ ~ 

n(x)  ~  c^  e  °  =  c^  e  e  ,  (3.50) 

which  is  a  Maxwell-Bol tzmann  distribution  in  the  energy. 

This  justifies  the  use  of  the  Maxwell-Boltzmann  distribution 
to  calculate  the  fluctuations  in  x  (3.29). 

To  study  the  time-dependent  solutions  of  (3.38)  we 

assume  that  the  variables  can  be  separated  and  hence  employ 

_  .  15,19 

the  ansatz 

n(x,t)=e  g(x)  ,  (3.51) 


iU.t) 
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which,  upon  substitution  into  (3.38),  gives 

Jr  V. 


1  9  rx+b  9g.  d  r  11, 

26u  9x  1 — x  9x  +  9x  a  +  26u~  g+  Ag  "  0 

u  o 


(3.52) 


Inserting 


g ( x)  =  (x+b) 


[26u  (b+a) (1+b) -1]  -26u  (b+a) (x+b)  6u  (1-x) 


o 


o 


h(x) 

(3.53) 


into  (3.49)  yields, 
26u  (b+a) 


_9_ 

9x 


(x+b) 


26u  (1-x) 
o 


-26u  (b+a) (x+b)  Su  (1-x)  . 

o  o  9  h  ( x) 

e  e 


-26u  (b+a) (x+b) +6u  (1-x) 

+  A  e  °  °  (x+b)  ° 


9x 


[26u  (b+a) (1+b) -1] 


h(.x) 


=  0 


(3.54) 


This  is  the  standard  form  of  a  Sturm-Liouville  equation. 
Hence  we  may  write  the  solution  to  (3.38)  as 


n (x, t)  =  e 


6u  [(1-x)  -2 (b+a) (x+b) ] 


(x+b) 


[26 u  (b+a)  (1+b) -1] 
o 


oo 


-A,  nt 

Tee  h  (x) 

L  n  n  n 

-n=0 


(3.55) 


where  A^  and  h^ (x)  satisfy  equation  (3.54)  plus  boundary 
conditions.  For  isothermal  desorption  the  appropriate 


boundary  conditions  are. 


n (x=l , t)  =  0  , 


and 


1  9n(x,t) 


1-x 


9x 


=  0 


(3.56) 


(3.57) 


x=l 
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Unfortunately  the  eigenvalue  problem  defined  by  (3.54), 
(3.56)  and  (3.57)  has  yet  to  be  solved.  However,  following 
Kreuzer  and  Teshima^^,  we  can  find  an  approximate  solution 
of  (3.38),  plus  boundary  conditions  (3.56)  and  (3.57),  by 
making  an  ansatz 

-t/t,  6u  (1-x) 2 

n(x,t)  =  Ae  a  (e  °  -  1)  ,  (3.58) 


which  assumes  a  single  desorption  time,  t^,  and  an  equili¬ 
brium  x  dependence.  We  also  ignore  terms  of  the  order  of 
the  thermal  fluctuations  in  (3.38)  and  get 

1  x  3n 


5n  _  _9_ 

3t  n  3x  l26u  1-x  9x 


+  xn 


(3.40) 


Substituting  (3.58)  into  (3.40)  and  cancelling  common 
factors  gives. 


5u  (1-x) 

( e  °  -  1 )  =  n 


-1  ,  o 


(3.59) 


Now  integrating  both  sides  from  zero  to  one  yields,  after 
rearranging  terms. 


6u  (1-x) 
e  °  dx  -  1 


Su 


0 


2g6u 

o 


72tt  WD 


-1 


2  6u 


m 


5u' 


(3.60) 


which  was  obtained  by  Kreuzer  and  Teshima.  The  dependence 
on  n,  which  has  already  been  obtained  by  Kramers^,  implies 
that  the  desorption  process  takes  longer  for  smaller  values 
of  the  friction.  Smaller  friction  results,  c.f.  (3.39), 
when  the  range  of  the  surface  potential,  ,  increases. 
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Although  (3.60)  was  previously  obtained  in  ref.  [10], 
the  above  calculation  did  serve  to  justify  some  of  the 
approximations  made  in  the  ansatz  (3.58).  In  particular, 
the  exponential  term  is  appropriate  because  all  terms  in 
the  general  eigenfunction  expansion  (3.55)  decay  exponen¬ 
tially  with  time.  For  our  analysis  we  also  argued  that  the 
quasi-equilibrium  assumption  implies  that  terms  of  the 
order  of  the  thermal  fluctuations  must  be  ignored. 

An  improved  approximation  for  the  desorption  time 
might  still  be  attainable  though.  Because  equation  (3.38) 
is  of  a  relatively  simple  form  we  hope  that  keeping  the 
thermal  fluctuations  may  give  a  more  accurate  result  for 
t^.  Alternatively,  the  lowest  eigenvalue  Xq  ~  might  be 
obtained  from  perturbation  theory. 


IV.  SUMMARY  AND  CONCLUSION 


A  detailed  comparison  of  the  low  friction  Kramers' 
equation,  as  derived  from  considerations  based  on  classical 
mechanics,  and  the  GKTT  equation  based  on  quantum  statis¬ 
tical  mechanics,  has  been  made.  This  was  achieved  by 
casting  both  equations  in  terms  of  scaled  action  variables, 
and  approximating  the  moments  of  the  GKTT  equation.  It  was 
found  that  if  terms  of  the  same  order  as  the  thermal  fluc¬ 
tuations  were  ignored  then  the  GKTT  equation  reduced  to 
Kramers'  equation.  In  contrast  to  the  derivation  of  the 
Kramers-Langevin  equation,  by  Caroli  et  al,  from  classical 

mechanics  in  the  limit  m  >>  M  ,  our  derivation  of  Kramers' 

s 

equation  holds  in  the  intermediate  temperature  regime. 


for  weakly  coupled  gas-solid  systems 


with  a  surface  potential  which  develops  a  large  number  of 
bound  states.  An  expression  for  the  friction  coefficient 
was  also  derived  and  was  found  to  be  identical,  up  to  a 
multiplicative  constant,  to  that  derived  by  Caroli  et  al. 
This  expression  shows  that  the  condition  of  weak  coupling 
is  indeed  equivalent  to  low  friction.  The  correspondence 
between  classical  and  quantum  theories  we  attribute  to  the 
large  number  of  one-phonon  events  necessary  for  desorption 
to  take  place  in  weakly  coupled  gas-solid  systems. 


Some  problems  on  the  subject  of  Fokker-Planck 


equations  for  desorption  remain.  The  derivation  of  the 
Fokker-Planck  equation  from  a  master  equation  has  yet  to 
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be  entirely  justified.  In  particular,  the  truncation  of 
the  Kramers-Moyal  expansion  for  a  system  with  finite 
boundaries  has  yet  to  be  tied  in  with  the  formalism  of 
van  Kampen's  large  parameter  expansion.  It  is  also 
desirable  to  solve  the  eigenfunction  and  eigenvalue 
problem  associated  with  the  GKTT  equation  so  that  a  more 
accurate  expression  of  the  desorption  time  may  be  found. 
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